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Abstract: The theoretical approach to describe the defocusing microscopy 
technique by U. Agero et al. [Phys. Rev. E 67, 051904 (2003)] assumes 
that the size of the objective lens aperture is infinite. This treatment 
gives that the intensity at the image plane depends on the laplacian of 
the phase introduced in the field by a pure phase object. In the present 
paper, we consider an arbitrary size for the aperture of the objective lens 
and we conclude that the intensity at the image plane depends also on 
the gradient of the phase introduced by the object and the phase itself. 
In this case, even an object that introduces only linear variations in the 
phase can be detected. Furthermore, we show that the contrast of the im- 
age of the phase object increases with the use of smaller objective apertures. 
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1. Introduction 

Understanding how a defocused lens system affects a diffracted image has been the theme of 
several investigations [1, 2, 3, 4, 5, 6, 7, 8, 9]. The experimental techniques and the developed 
theoretical background has been used as a tool for obtaining spatial and volumetric information 
of a microscopic phase object [5, 6, 7, 8, 9, 10, 11, 12]. Throughout the last 60 years many 
other microscopy techniques were developed for this purpose, such as the well-known Zernike 
phase-contrast method [13], the digital speckle photography [14], diffraction phase microscopy 
[15], diffraction tomography [16], dark field microscopy [17], differential interference contrast 
microscopy [17], ghost images [18] and techniques which use the chromatic aberration [19]. 
However, in comparison with these techniques, the defocusing microscopy technique (DM) is 
the simplest to implement, since we just need displace in the objective lens of the microscope 
from the image focusing distance. 

Despite of its simplicity, DM is a very powerful technique, which has been applied in the 
research of biological systems. It was used to observe the shape, speed and density of coherent 
propagating structures in macrophages (cell of the innate immune system), and measure the 
time of a single phagocytosis event [9]. Moreover, the shape, size, refractive index, bending 
modulus, and cytoplasm viscosity of red blood cells was obtained by applying DM to these 
cells [20]. DM was also used to do the tomography of fluctuating cell surfaces [11]. The cur- 
vature of the cell membrane is directly related to the laplacian of the phase introduced by it 
in the propagating field. On the other hand, the laplacian of the phase can be obtained from 
the contrast of the image observed when DM is used [9]. Thus, DM allows for obtaining the 
curvature of the surface of a cell directly, avoiding a subsequent calculation of the curvature. 

Initially, the term "Central illumination microscopy" was used to refer to DM [2], where it 
was applied for determining the refractive index of transparent materials (refractometry). In 
this case, one usually observes an interference pattern which has been called "Becke line" [21]. 
This pattern is symmetric when the phase edge varies abruptly and asymmetric when the phase 
varies slowly. 

The theoretical approach for DM of Ref. [9] assumes that the size of the objective lens 



apertures is infinite. This approach predicts that the intensity at the image plane of a pure phase 
object depends only on the laplacian of the phase being modulated in the field. So, any linear 
phase change introduced by the object should not be detected by DM. In this paper we study the 
use of arbitrary apertures sizes for the objective lens and we show that by using small lenses, 
the intensity at the image plane depends also on the gradient of the phase and the phase itself. In 
this case, even an object that introduces only linear phase changes in the field can be detected. 
This improvement may lead to new applications for DM, such as the observation of "hidden" 
structures in cells. 

Moreover, we study the dependence of the image contrast with the size of the objective 
lens aperture. We perform an experimental test of our theory and show that the image contrast 
increases when the entrance pupil of a microscope system is reduced. It is important to observe 
that for DM, to the best of our knowledge, there is no investigation about the role of the finite 
size of the lens aperture for the case of a coherent illumination. Only in the case of incoherent 
and partial coherent illuminations there have been studies showing that smaller lenses increases 
the contrast of the image due to the coherence introduced by spatial filtering [2, 6]. 

The rest of the article is structured as follows. In the following section, we make a short 
review of lens systems for imaging formation. In section 3, we present the main theoretical 
results of this work, which is the extension of DM technique for the case of a finite objective 
lens aperture. In section 4, we show a DM experiment in which it is possible to observe the 
effect of the finite size of the aperture in the contrast of the image of the phase object. It is 
exposed the concluding remarks of this manuscript in section 5. 

2. Lens systems for image formation 

2.1. The simplest lens system 

The simplest imaging system uses a single lens, as shown in Fig. 1(a). S is a source of plane 
wave field. LI is a lens, O is the object plane, and I is the image plane. In this case, the ex- 
pression for the intensity of light at the position p =x i + y j in the image plane, supposing a 
coherent illumination, is the square modulus of the following electric field [22] 

E(p)=Ax Jd^Eo^)T(^ + ^j, (1) 

where Z\ is the distance from the object plane to the lens and Z2 is the distance from the lens 
to the image plane, k is the magnitude of the wave vector. £0 (£) is the light electric field at 

the position | = ^ x ~t + E,y~~j immediately after the object plane. T + is the Fourier 
transform of the transmission function of the lens aperture. The term A before the integral in 
Eq. (1) contains only scale factors and phase components that do not affect the form of the 
image intensity. All variables which appear in boldface in all equations of this paper mean 
vectors in the transverse plane having x and y components. The Eq. (1) is valid whether one of 
the following two conditions are satisfied: (i) the object is illuminated by a spherical wave that 
is converging towards the point where the optical axis pierces the lens [22] or (ii) the size of 
object is no greater than about 1/4 the size of the lens aperture [22, 23]. 

2.2. A simple microscope system 

A slightly more complex imaging system discussed in [9] is the microscope constituted by two 
lenses, as can be seen in Fig. 1(b). Considering a finite aperture for the objective lens and an 
infinite aperture for the eyepiece lens, the expression for the intensity of light in a position p at 
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Fig. 1. Schematic draw of three image systems. In these figures, S is a plane wave source 
that illuminates the object, represents the object plane, I represents the image plane, and 
f\ and fo are the focal lengths of the lenses Ll and L2, respectively, (a) Simplest system to 
observe the image of an ordinary object. Object and image planes are placed at distances 
Zj and Z2, before and after the lens Ll, respectively, (b) Microscope system. This system 
is used to generate the image of ordinary objects. The distance between lenses Ll and L2 
is equal the sum of their respective focal lengths, (c) Defocusing microscopy scheme. Af 
is the length of the defocusing. This scheme is based on the one considered in Ref. [9]. 



the image plane, supposing a coherent illumination, is given by (see appendix) 

E(p)=Bx j d $E {$)T i (^j- + k j^, (2) 

where f\ and fi are the focal lengths of the lenses considered in the setup. T\ + j£\ is 
the Fourier transform of the transmission function of the objective lens aperture. B is a term 
that includes only scale factors and phases. Since we are interested in the intensity detected 
at the image plane, from now on we will no longer write this term in order to simplify the 
mathematical manipulation. Note that to form an image, it is assumed for Eq. (2) the same 
conditions considered for Eq. (1) mentioned in the last section [22, 23]. 

2.3. Ideal imaging systems 

Better image resolutions can be obtained by lenses with large apertures. By assuming infinite 
apertures in both Eqs. (1) and (2), their Fourier transforms become delta functions which filter 
the integrals. In the first case [Eq. (1)] we obtain 

E(P) = E o(~p\ (3) 



In the second case [Eq. (2)] we obtain 



E(P) = M-J/)- W 
In both cases, the minus sign indicates that the image is inverted when compared to the object. 
2.4. Image formation of pure phase objects without special methods 

For a pure phase object and uniform illumination, we have Eq (§) = e l ™\ which results that 
the square modulus of both Eqs. (3) and (4) are uniform for infinite aperture lenses, since this 
last condition means that T and T\ are Dirac delta functions. In the case of an amplitude object, 
a decrease of the size of the lens aperture implies in a decrease of the resolution of the image 
[24, 25]. So, one could expect that the use of lenses with smaller apertures would make the 
methods described in sections 3 and 4 even more useless for the image formation of pure phase 
objects. Fortunately, this is an incorrect assumption (see Ref. [21]). In this section we discuss 
an example in order to show this. 

Let us consider a pure phase object that introduces a n of discontinuity in the phase of the 
incident field. This discontinuity is located at the position 1 mm from the optical axis of prop- 
agation of the field at the plane of the object, as can be seen in Fig. 2(a). Figures 2(b)-(g) 
show the intensity of the light at the image plane. These intensities are calculated by using Eq. 
(2) and considering the widths of the cylindrical lenses equal to 2L = 300 mm, 2L = 200 mm, 
2L = 100 mm, 2L = 10 mm, 2L = 5 mm and 2L = 1 mm, respectively. For simplicity, we con- 
sider an one-dimensional case. We suppose the use of a bi-concave cylindrical lens, that pro- 
vides magnification only in one dimension. This implies that the one-dimensional Fourier trans- 
form in Eq. (2) is given by sine + ■ 

Note that the fringes around the position —5 mm are owing to the phase discontinuity in 
the position 1 mm at the object plane. These positions are related by the magnification M = 
—fllf\ = —5. It is important to note that due to the small size of the objective it is possible 
to detect the discontinuity by means of the usual method of image formation, without any 
defocusing or introduction of some other contrast technique. That is a result that can not be 
obtained if infinite lenses are considered. In Fig. 2 we can also note that the diffraction bands 
at the image plane become narrower when the size of the objective lens increases, which is in 
agreement with Ref. [2] . 

Figure 2(h) shows the maximum variation of the intensity in the image plane with respect to 
different widths of the cylindrical lenses. This variation is defined as AI max = I max — I m j n and it 
is directly related to the contrast of the image [9]. Figure 2(h) shows that objectives with small 
apertures (i.e. smaller values of L) generates images with higher contrasts. This behavior means 
that a restriction in the linear dimension of the lens aperture improves the capability to detect 
phase objects. 

3. Defocusing microscopy 

3.1. Defocusing microscopy with an objective lens of arbitrary size aperture 

The DM consists of making a small displacement at the objective longitudinal position, as 
shown in Fig. 1(c). This displacement allows the detection and reconstruction of phase objects. 
The treatment of Ref. [9] assumes that the lenses apertures are infinite. This approach gives the 
following expression for the intensity at the image plane 
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Fig. 2. (a) Phase introduced by a pure phase object. It has a discontinuity at the position 
1 mm. (b)-(g) Intensity of the images produced by cylindrical lenses with different widths. 
The intensities, which were obtained assuming an ordinary microscope system (i.e. without 
defocusing), were calculated by using Eq. (2). The parameters used in the simulations are 
f\ = 100 mm, fo = 500 mm, and 826 nm of wavelength. The widths were 2L = 300 mm, 
2L = 200 mm, 2L= 100 mm, 2L = 10 mm, 2L = 5 mm and 2L = 1 mm in figures (b), (c), 
(d), (e), (f) and (g), respectively, (h) The maximum variation of the intensity {M max = 
Imax — Imin) in figures (b)-(g) as a function of the width 2L. 



where p' = — j^p and the subscript p' indicates that the laplacian is calculated with respect to 
this variable. The first term of Eq. (5) is associated to the backlight, which is always present in 
the image formation of phase objects. The parameter Af in the second term is the amount of 
defocusing [see Fig. 1(c)]. It is important to observe that any linear change in the phase is not 
taken into account by Eq. (5), owing to the fact that the laplacian of a linear function is zero. 

3.2. Defocusing microscopy with an arbitrary size for the objective lens aperture 

In this section we will use the Fresnel propagator as a tool in order to study the case of DM 
with finite size for the objective lens aperture. Consider the microscope system depicted in Fig. 



1(c). The electric field at the image plane is given by 

E(p) = JdtEoffiJdafaWle^effi^n 

dj5|A 2 (j5)|e^" cW e ~^h , (6) 

where £o (§) is the electric field immediately after the object plane. |Ai(ff)| and |A2(/5)| are 
the modulus of the transmission function of the objective and eyepiece, respectively. f\ and fi 
are the focal lengths of the objective and eyepiece lenses, respectively. The integration in the 
variables |, a and /3 are due to the propagation from the object to the objective, objective to 
the eyepiece and eyepiece to the image, respectively. Eq. (6) can be written as 

ikp 2 r r 



where 



E{p)^e~ 1 ^ / dl;Eo(£) / c/a|A 1 (a)|e 2 W- A - f )e^7Te 2 (''+ A /)/j(a,p), (7) 



r -ikp 2 ik[p 2 -2pa] ik\p 2 -ipp\ 

h(a,p)= dfieVTe e J h . 



(8) 

In this last step we assumed that the eyepiece is sufficiently large so that |A2(/3)| ~ 1. The 
exponentials depending on /}" cancel out, and the expression (8) simplifies to 

-ik <„ , (rf+A/) l2 

/i(a,p)=CxeW) l + ^ PJ , (9) 

where C is a term that represents only scale factors that does not depend upon the variables 
CC and p and will be dropped out. Inserting Eq. (9) into Eq. (7), expanding all quadratic ex- 
ponentials, and performing the algebraic sum of quadratic terms in the variable Ot, we obtain 
that 

,k\ i _f ll frW] 2 ik , 



^a|A 1 (a)|^(^-7r) 0t2 e -''K7T i A7 + 7|)- tt . (10) 



Note that Eq. (10) is valid even for large values of Af. Now let us suppose that Af is small, 
so that we can write 7^57 ~ J[ + ( m i s expression is obtained by Taylor expansion around 
Af = 0). Thus, the first exponential expression in the second integral of Eq. (10) becomes 

Unless a phase factor (which is irrelevant to the calculation of the intensity), the eletric field is 
given by 

E(p) = J d^Eoi^e 2 ^^ J dalA^l^l + '^^y-^T^+Tlh. (i 2 ) 
By using the identity 



it is possible to show that 
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where is the Fourier transform of |A i (ce) | . Now let us define the function E' (p) as 
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Note that \Eq(p)\ 2 give us the image of the phase object generated by a microscope system 
whose magnification is M = j^^j ~ It is valid for a coherent illumination or for a partial 
coherent illumination when the object is smaller than the transverse coherence area of the field, 
as pointed out in Ref. [9]. From now on, it will be assumed that one of these two physical 
conditions is respected. By using the above definition in Eq. (14), we have 



E{p)=E' {p) -^V 2 ^(p). 



(16) 



This expression contains the laplacian of the Fourier transform of the modulus of the objec- 
tive lens transmission function, the only part dependent on p . Let us try to get an expression 
containing the laplacian of the object function [as in Eq. (5)]. By remembering that 
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(18) 



This expression contains the laplacian of the function of the object, however £0 is calculated 
at the coordinate {f\ — A/)(| — ■£■) instead of the coordinate v. On the other hand, the Fourier 

transform which is usually calculated at the coordinate {j- + £-) [see Eq. (2)], here is evaluated 
at v. In order to write the above expression in a simpler form, let us make the change of variable 
given by 



C/i-4fl(T -£) = «■ 

k 12 



(19) 



kt 1 kp 



This change implies that v = j^t^j + , dv = j^^jd^ and V 2 = V 2 , 
the subscripts indicate the variables for which the laplacian is calculated. Using this change of 
variable and neglecting the terms that contain powers of A/ larger than or equal 2 (as it was 
done in [9]), we obtain that 
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Eq. (20) is valid for small defocusing and any general class of objects. Let us now suppose a 
pure phase object illuminated by a plane wave with amplitude Equ such that the transmitted 
complex field amplitude is given by 

E (Z)=E 0i e i *®. (21) 
We also assume that the size of the object is smaller than about 1 /4 of the size of the lens 



aperture, so that we can neglect the quadratic phase factors e 2 (/i- A /) according to [22, 23]. In 
this way, it follows that the electric field is given by 



E(p) = J d$E 0i 
By using again the identity 
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we obtain 
E{p) = 



V 2 (e^)=e^ [-(V0) 2 + /V 2 </>] , 
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This is the final expression for the electric field. The intensity at the image plane, I(p) 
\E(p)\ 2 , is given by 



I(p) = fd^Eoie^ 
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Note that it is not necessary to take the conjugate of the function T\ because it is a real quantity. 
By multiplying the three terms in each integral for their conjugates we obtain nine terms. By 
neglecting all quadratic powers of Af and considering EoiE^ a constant factor we obtain the 
final expression for intensity 
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This is a real function. The approximation 
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above implies that we may replace f\ — Af by f in the denominator of the argument of T\, 
which gives 
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When the size of the lens aperture tends to infinite we have T\ y + ^ 
where a — £, , 77 .In this case, the delta functions filter the integrals, so that 
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Fig. 3. (a) Object complex transmission function, (b) The intensity at the image plane when 
2L = 10 mm and Af = 0. (c) The intensity at the image plane of the same phase object, but 
now considering Af = 2 mm and 2L = 200 mm. (d) The intensity at the image plane for 
2L = 10 mm, keeping Af = 2 mm (See the main text for details). 



which reproduces Eq. (5) as a particular case. Also, p' = — y-p and the subscript p' indicates 
that the laplacian is calculated with respect to this variable. 

Equation (27) is the main theoretical result of this paper. It is a generalization of the previous 
treatment given for DM in Ref. [9]. Note that the gradient and the phase itself are now present 
in the expression for the intensity in the image plane. So, an object that introduces only linear 
changes in the phase of the field can be detected by DM with the use of a small objective 
aperture. This is not possible when an infinite objective aperture is considered. 

Figures 3(b)-(d) show the intensities at the image plane for a pure phase object that introduces 
a linear phase change into the field, as shown in Fig. 3(a). In this case, Eq. (5) predicts only 
an uniform backlight. These curves were calculated via Eq. (27) with the parameters f\ = 
100 mm, f2 = 500 mm, and 826 mm of wavelength. Figure 3(b) shows the intensity at the image 
plane when 2L = 10mm and Af = 0. The oscillations around the points p = ±0.5 mm at the 
image plane are related to the points where there are abrupt phase changes at the object plane, 
considering the magnification of the imaging system, M = — /2//1 = —5. In Fig. 3(c), we have 
the intensity at the image plane of the same phase object, but now considering Af = 2 mm. It 
was also assumed an objective lens aperture of 2L = 200 mm. The idea of using a large objective 
aperture is to simulate an infinite lens. Note that the values of f\ and Af are compatible with 
the approximations j^^j ~ j- { + in the quadratic phase factors and /1 — Af « f\ in the 

denominator of the argument of T\ . The contrast of the diffracted images tends to zero when 
the size 2L of the objective lens tends to infinite. This shows that a small objective aperture can 
be more effective than defocusing the system in some situations. Finally, Fig. 3(d) shows the 
intensity for a lens aperture of 2L = 10 mm, keeping Af = 2 mm. For this case, the intensity 
has more contrast in comparison with the last ones due to the joint action of the finite size of 
the objective lens aperture and the defocus of the objective lens. 

4. Experiment 

We performed an experiment to verify some predictions of our theoretical generalization of 
DM, presented in the last section. More specifically, we experimentally study the dependence 
of the image contrast with the aperture of the objective lens of a microscope system. The ex- 
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Fig. 4. Schematic draw of the experimental setup. L\ is a 25 mm focal length lens. The focal 
length of L2 is 200 mm. The focal lengths of L3 and L4 are both 300 mm. L3 is mounted on 
a motorized linear stage. PH is a 100 flm diameter pinhole. P, (i = 1 , 2) are polarizers. SLM 
is the LCD of the phase spatial light modulator. M, (; = 1,2, 3) are mirrors. Diode laser is a 
633 nm laser. CMOS is a camera. All lenses have 1" in diameter. 

perimental setup is showed in Fig. 4. A single mode diode laser operating in TEMoo mode with 
~ 0.5 mm FWHM and at 633 nm wavelength was used as the source of illumination. This beam 
crosses a telescope lens system with a spatial filter to improve its quality and the beam spot size. 
After the telescope, the beam is well collimated and magnified by a factor 8. We used lenses 
LI and L2 (1" in diameter) that have 25 mm and 200 mm of focal lengths, respectively. We also 
used a pinhole (PH) of 100 jJ.m of diameter that was placed between the lenses. This allowed us 
to obtain the spatial beam profile shown in Fig. 5. It shows that for illuminated objects smaller 
than 1 mm, it is reasonable to assume that the incident field is a plane wave. The beam is sent 
to a reflective spatial light modulator (SLM), which is composed of two polarizers (PI and P2) 
and a twisted nematic liquid crystal display (LCD), model Holoeye 1080P. The SLM is adjusted 
for phase modulation only [26]. 




"1000 -3000 -2000 -1000 1000 2000 3000 4000 

Position (urn) - horizontal direction 



Fig. 5. Beam profile (at the horizontal direction) of the laser after crossing the telescope. 
The data was obtained with a scanning slit optical beam profiler. 

Figure 6(a) shows the curve of the phase modulation versus the LCD gray level for our SLM 
at 633 nm wavelength. Using the SLM at the phase modulation configuration, we tested if there 
was still any amplitude modulation in the incident field. The amplitude as a function of the LCD 
gray level is shown in Fig. 6(b). This curve ensures that the equipment is modulating the phase 
of the field with a fluctuation of less than 10% in the amplitude modulation. After ensuring the 
capability to generate a reasonable pure phase object, we generate at the SLM an object that 
introduces a phase discontinuity at the incident field equal to n [See Fig. 6(c)], as shown in Fig. 



6(d). For simplicity, we considered a one-dimensional phase object. It simulates the effect of 
the edge of a transparent object over a plane wave field, supposing that this field has crossed the 
object in a perpendicular direction to its plane face [2, 6] . In fact, since a pixel of this equipment 
has 8 jJ.m of linear dimension, and also because the SLM is based on a reflective LCD, we can 
guarantee that the phase discontinuity occurs in a very small zone over the horizontal direction. 
To perform a n phase change, we have chosen for the LCD gray level at the left side of the 
discontinuity and 52 for the LCD gray level at the right side. The object was generated at the 
SLM such that the discontinuity was located at the center of the incident pump beam. 




Fig. 6. (a) Calibration curve of the SLM for phase modulation at 633 nm wavelength, (b) 
Amplitude modulation versus the LCD gray level, (c) Modulation used and (d) the corre- 
sponding generated object at the SLM. 



After being reflected by the LCD, the light follows to a microscope system, which is prepared 
with the lenses L3 (objective) and L4 (eyepiece), as can be seen in Fig. 4. They have 1" of 
diameter and equal focal lengths of 300 mm. In this system, the LCD is placed at the focal 
distance of L3, L4 is placed at a distance of 600 mm from L3, and the camera is positioned at 
300 mm after L4 at the image plane. The rectangular variable aperture (VA) is placed right in 
front of the lens L3 in order to control the dimension of the entrance pupil of the system (see Fig. 
4). It has a fixed length in the vertical direction and variable length in the horizontal direction. 
Its horizontal length is controllable by a screw with a 0.005" of precision. The VA together 
with L3 constitute an effective objective lens with a controllable horizontal linear aperture. To 
control the defocus, the system VA+L3 is placed over a motorized stage that allows for the 
control of its longitudinal position with a precision of 100 fim. The camera used to obtain the 
images is a CMOS, model EO-1312M 1/2" Monochrome USB Lite Edition, from Edmond 
Optics. 

Figure 7 shows the obtained experimental results. The interference patterns were recorded 
with Af = 0, in order to verify the effects associated with the variation of the size of the entrance 
pupil only. Figures 7(a), 7(b), and 7(c) show the intensities of the noise-subtracted images, 
while considering VA equal to 0.0600", 0.130", and 0.190", respectively. Given an aperture, 
two images were measured. They correspond to the case where the object was being or not 
generated in the SLM. When SLM generates the object, the image has contributions from the 
noisy background and from the phase modulation performed by the device. This means that 
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Fig. 7. (a), (b), and (c) show the intensities of the noise-subtracted images versus horizontal 
position in the plane of the camera for apertures equal to 0.0600", 0.130", and 0.190", 
respectively, (d) shows AI mm versus the linear size of the objective lens aperture. 



it is necessary to subtract the noise from the image to proper take into account the effect of 
the entrance pupil in the Becke line. To do this, we recorded the contribution to the image 
corresponding to the background noise. This is done by generating an uniform global phase 
object at the SLM, such that it works as a conventional mirror. This procedure reveals all non- 
desired modulations at the image plane. Then, we subtract this noisy image from the signal 
recorded before. An example of a noise-subtracted image can be seen in Fig. 8. Note that the 
subtracted intensity (AT) can be negative since the intensity of the noise in a pixel can be greater 
than the diffracted image when the SLM is generating the object. 

In order to evaluate A/, we performed a statistical analysis of each noise-subtracted image 
obtained. Each image was treated as a matrix, where each element was associated with a phys- 
ical pixel of the camera. The corresponding AI for a certain horizontal position at the image 
plane was calculated by taking the mean value of the intensities of each pixel in a column. Each 
column was associated with a certain horizontal position. As one can see from Figs. 7(a), 7(b) 
and 7(c), the contrast increases when the linear horizontal dimension of VA is reduced. Fig. 7(d) 
shows the curve for AI max as a function of the apertures measured. It is possible to associate 
Almax directly with the contrast of the image observed [9], where larger AI max implies higher 
contrast. 

5. Conclusion 

We have studied theoretically and experimentally the defocusing microscopy technique by tak- 
ing into account the finite size of the objective lens aperture. We concluded that by using smaller 
objective lens aperture one can detect linear phase changes caused by phase objects, that could 
not be detected when infinite (very large) lens are considered (used). This improvement may 
lead to new applications for DM, such as the observation of "hidden" structures in cells. 

We also observed that small lens apertures generates diffracted images with higher contrast 
than lenses with infinite aperture diameters.. This type of manipulation can be combined with 
defocusing to improve the efficiency of the DM technique. It is important to note that DM is a 




Fig. 8. The noise-subtracted image for the horizontal aperture of VA equal to 0.100". Each 
pixel of the CMOS has 5.2 /lm. 



low resource consuming technique that can easily be applied in any microscope. Thus, the work 
presented here will be relevant in the study of any phase object in different areas of science. 

A. Appendix: Microscopic System 

Using the Fresnel propagators [22], we found that the electric field at the image plane for the 
apparatus shown in figure 1(b) is given by 

E(p) = Jd$E (Z) J da |Ai(a)|«~3re~*r- 

r -ikfp tt|P-a| 2 ifr|p-p| 2 

d^\A 1 {^)\e'^~e^h+f2)e~^~ , (29) 

where |A i (ot) and |A2(j3)| are the modulus of transmission function of objective and ocular 
lens respectively. The integrations in Ct and j3 implement the propagations from the ob- 
ject to the objective, from the objective to the ocular and from the ocular to the image plane 
respectively. The Eq. (29) can be written as 



(30) 



where 



r r ik\p-a\- -agj -il:pp 

g(£,p) = J daj dfi\A 1 {a)\\A 2 {fi)\e^^e~fi~e-S~. (31) 

Let us define 

k% kp k 

J =Xl ' 7T* 2 ' 2(KTW)= y - (32) 

In this way we have 

g(x u x 2 ) = FF{\A 1 (a)\\A 2 (P)\ e '^-«l 2 }, (33) 

where the symbol FF mean the double Fourier transform with respect to the variable par (ot, /3 ). 
Using the convolution theorem [22] we have 

g(*l,* 2 )=FF{|Ai(a)||A 2 (0)|}* FF{ e '-^- a l 2 }, (34) 

where the symbol * means the convolution. The double Fourier transform assumes the fac- 
torable form 

FF{\A l (a)\\A 2 (m = F{\A 1 (a)\}F{\A 2 (p)\} = T l (x 1 )T 2 (x 2 ). (35) 
Let us define the function h(xi , x 2 ) as 

/ 1 (jc 1 ,jc 2 ) = FF{e ,v ^- a l 2 } (36) 

so that 

g(xi,x 2 ) = Ti(xi)T 2 (x 2 ) * h(x u x 2 ). (37) 
After some calculations it is possible to show that 

h(x u x 2 ) = 8{x x +x 2 )e^ £L $ L . (38) 

In order to demonstrate the above equation we suggest the following change of variables: /} + 
(X = 2u and j3 — (X = 2v. These changes will eliminate the cross terms so that the expression 
will become the product of two independents integrals. Inserting Eq. (38) in Eq. (37) we have 

g(x h x 2 ) = T 1 (x 1 )T 2 (x 2 )*8(x 1 +x 2 )e-'^ L 1 (39) 

that is, 

r r - i4 -A i 2 

g(x u x 2 ) = / dx\ / dx' 2 Ti(xi -x\)T 2 {x 2 -x' 2 )8{x\ +x' 2 )e '~ ^~ . (40) 
Solving the above integral and going back to the original variables we obtain 

*(4,p) = JdvT l (^-- V ^jT 2 (^+ V ^e^ l2lv2 . (41) 



This is the optical transfer function (OTF) of the microscope system, which takes into account 
the finite size of both objective and ocular apertures. Assuming an infinite ocular aperture, so 
that To + v^J = 8 + v^j we conclude that 

gg,p)= T Aj- + j-y 2 h . (42) 
Inserting the Eq. (42) in Eq. (30), we find 

E{p)= const, e^-e 2f 2 / d$E (€)e^7i ^- + . (43) 
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According [22], the quadratic phase factor e"W in Eq. (43) causes "an unacceptably large 
image blur". In order to eliminate it, we consider one of the following two restrictions. First, 
we can assume that the object is illuminated by a spherical wave that is converging towards the 

-i»g 2 

point where the optical axis pierces the lens, that is, the incident field contains the term e 2/ i 
which cancels the quadratic exponential [22]. Second, we can assume that the size of object is 
no greater than about 1/4 the size of the lens aperture [22, 23]. After assuming one of these 
restrictions, we have finally the expression 

, kp 2 -■*(/! +/ 2 ) 2 „ fkE kO\ 

E{p) = const, etc 2 /| / d%E (€)7i I ^ + ) . (44) 



/i h 

The phase factor outside the integral vanishes only when the size of both (ocular and objective) 
lenses tends to infinity. 
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